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Analytical and numerical studies on many-body stochastic processes with multiplicative 
interactions are reviewed. The method of moment relations is used to investigate effects 
of asymmetry and randomness in interactions. Probability distribution functions of the 
processes generally have similarity solutions with power-law tails. Growth rates of the system 
and power-law exponents of the tails are determined via transcendental equations. Good 
agreement is achieved between analytical calculations and Monte Carlo simulations. 

§1. Introduction 

For several decates, power-law distributions have been one of the central issues in 
statistical physics and have been observed not only in physical, but also in biological, 
social, and economic phenomena. In order to explain this spontaneous emergence 
of power laws, several concepts, such as self-organized criticality,^-*"^^ scale-free net- 
work,'^) and so on, have been proposed. Some power laws have been successfully 
explained by the concepts, but others are not sufficiently. A lot of researchers have 
strenuously engaged in finding underlying physics of the power-law distributions. 

In recent years, inelastic Maxwell models dealing with the system of inelasti- 
cally colliding particles have been investigated and have revealed that the veloc- 
ity distribution function of particles has a power-law tail.^^'^^ Inelastic Maxwell 
model(IMM) is defined on A^— particle system. Each particle has a positive quantity 
Xi (i = 1, . . . , N) and at each timestep, two of them are selected randomly and expe- 
rience a binary interaction for transforming the two quantities of particles Xi , Xj into 
= exi + (1 — e)xj, x'j = (1 — e)xj -|- exj where < e < 1 is a restitution coefficient. 
Similar models such as greedy mutiplicative exchange(GME)^^ whose interaction is 
described by = (1 — a)xi, x'j = Xj + axi {xj > x-i) where < a < 1 is a interaction 
parameter and symmetrically multiplicative interaction(SMI)^-* whose interaction is 
x[ = axi + Pxj, x'j = f3xi + axj where a,(3 > have also been studied. Character- 
istic features common to these models are that they have similarity solutions scaled 
by a growth rate of the system and power-law tails in the probability distribution 
function(PDF) of the quantities. In general, the power-law tails of these models 
appear if their interaction parameters are within a finite region rather than a critical 
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point. This means that fine-tuning of the parameters is not needed for observing the 
power laws, while it is necessary in usual critical phenomena. In addition, it is also 
found that the power-law exponents, which are able to be calculated analytically, 
vary continuously with the values of interaction parameters. 

In this paper, we collectively call these models many-body stochastic processes 
with multiplicative interactions and consider two extended models: asymmetrically 
multiplicative interaction and effects of randomness in interaction. The contents are 
organized as follows. Firstly, key points of theoretical procedures of analyzing the 
processes are briefly introduced. Secondly, a model of asymmetrically multiplica- 
tive interaction is considered. Thirdly, a model which takes into account effects of 
randomness in interaction is considered. Finally, we conclude with a summary and 
application of our models. 

§2. Key points of the theory of many-body stochastic processes with 

multiplicative interactions 

Let us introduce general procedures to investigate the tails of PDFs of the pro- 
cesses analytically. The method of Fourier transform has been often adopted for this 
issue.^^'^'*'^^"^'^^' However, we would like to consider it here with using a method 
of moment relations which is easier to understand it intuitively. We will apply this 
method to our models in the next two sections. In order to introduce the procedures, 
let us take SMI whose interaction rule is x' = ax + (3y, y' = fix + ay as an example. 
In the limit the number of particles ^ oo, a master equation of the system is 
descirbed as follows. 

O /'OO /'OO 

-f(z,t) = J^ dxj^ dy f{x,t) f{y,t) [5{z-{ax + (3y))-5{z-x)]. (2-1) 

The goal of our calculations is to find a similarity solution of this equation. Thus, 
the following scaling relations are assumed. 

i = ze-^\ f{z,t)=e-''m, (2-2) 

where 7 which is obtained explicitly later is a growth rate of the system. Substituting 
Eq. H2-2() into Eq. H2-l() and some calculations with using Lemma 2 in Bobylev et 
al.^s paper, the next inequalities of moments are obtained. 



Ekp 1 
k=l 



P7 - {aP + PP -1) 



Eti (l) {a^P^-'' + aP-^l3^),,ul^p^u 

< f^P < ^ ( 1^ ' (2-3) 

— [aP + pP — 1) 

where Hp is the p—th order moment of l^(^) and kp denotes the integer part of 
{p + l)/2 for p > 1. Note that when p is an odd integer, the first inequality becomes 
an equality. This determines the growth rate j = a + P — 1 at p = 1. The moment 
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inequalities (|2-3|) indicate that if a transcendental equation = + (3^ — 1 has 
a solution except a trivial one p = 1, let us say p = s, then the moments fip of 
(1 <) p < s are bounded and those of p > s (> 1) diverge. This immediately lead 
to a power-law tail with the exponent s in PDF l^(^) ~ 1/.^^^* where ^ ^ 1. When 
< s < 1, this formalism based on Eq. (|2-3p breaks down. Instead, an inequality 
P'j < + — 1 can be derived in this case. This indicates that a power-law tail 
appears if the system is scaled by the minimum growth rate jmin under the condition 
that s'jmin = + (3^ — 1 has a solution. Thus, the exponent s is determined by 
{q^ lna + (3^ ln/3}s = a' + f3' - I. 

Consequently, the power-law exponent s is determined by transcendental equa- 
tions {a+P-l)s = a'^+P'-l {s > 1) and {aMn a-F/3Mn/3}s = a'+f3'-l (0 < s < 1) 
and the growth rate becomes 7 = a-|-/? — 1 (s > 1) and 7 = ^min 

(0 < s < 1). 

This denotes that 7 and s are continuous functions of interaction parameters a and 
p. These are the key points of the theory of many-body stochastic processes with 
multiplicative interactions. 

§3. Effects of asymmetry in interaction 

In this section, we investigate many-body stochastic processes with asymmetri- 
cally multiplicative interactions ( AMI). ^'^^ A interaction rule of the processes is ex- 
pressed as x' = c(l — a)x + cby, y' = dax + d{l — b)y {x > y) where < a, 6 < 1 
and c, d > are interaction parameters representing exchange rates and amplifica- 
tion rates of larger and smaller quantities, respectively. This rule distinguishes one 
interacting particle from the other by the magnitude of thier quantities x, y. This 
difference of interaction rule between two interacting particles is what we call asym- 
metry. Note that this model includes the three models IMM(c = d = 1, a = 5), 
GMI(c = d = 1, a = 0), and SMI(c = d, a = b) aforesaid as special cases. We have 
performed Monte Carlo simulations of this model and have found that PDF obeys 
a power law at the tail as shown in Fig. ^i^a). So, we examine this power-law tail 
in accordance with the procedures mentioned in Sec. [2 A master equation of this 
system is described by 

+ f{z, t) = dy 1" dxf{x, t)fiy, t) 

X [6{z - (c(l - a)x + cby)) + 6{z - {dax + d{l - b)y))] (3-1) 

The same scaling relations H2-2() is also used in this case. Substituing Eq. ()2-2() into 
Eq. ()3-l() and some calculations lead to the following moment inequalities. 

E ( ^ ) [ - a)}'{cbr-' + {c(l - a)r-''{cb}' + {da}'^{d{l - 6)}^"^ 
k=i ^ ^ 

/•oo roo 

+{dar-Hd(.i - b)}' ] X / d^j d6(ei'er' + ei-'e2M^iM(2) 

Jo J £,2 

/"OO f'OO 

< (P7+iK-({c(i-a)r + {dan / di2 deieMeiM6) 

Jo J6 
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POO poo 

-{{chY + {d{i - h)Y) / d6 / 

- ( ) [ •^^(^ " + - o)Y-^{chY + {daY{d{l - 

fc=i ^ ^ 

Jo J^2 

Here, we assume that a scaled PDF obeys a power-law tail l^(^) ~ l/^^^'*, which is 
supported by the results of Monte Carlo simulations. It can be estimated that in 
the limit p ^ s, two integral terms in the above inequalities ()3-'2p satisfy the next 
relations. 

/»CXD poo poo poo 

/ dC2 / d^ieimiMC2) » / dC2 / rf^ef-f (6)1^(6), (3-3) 

Jo Ji2 Jo J(,2 

poo poo 

/ d6 / d^iei^{^i)^{^2) ^fip^oo. (3-4) 

Jo J^2 

This immediately follows that the growth rate of the system 7 and power-law expo- 
nent s of PDF are determined by the following transcendental equation. 

S7 = {c(l - a)y + {day - 1, (3-5) 
7 = ^ [c(l -a + b) + d{l + a-b)] + ^{c-d){l-a-b) A-1, (3-6) 

1 /"OO /'OO 

A = — di2 dCi (6 -6)1^(6)1^(6). (3-7) 

/^l Jo Ji2 

Consequently, we can successfully derive a transcendental equation (|3-5p which de- 
termines the power-law exponent s for s > 1. Note that the value of integral A in 
Eq. 1)3 -6^ is unable to be calculated analytically, thus a numerical integration must 
be performed to solve the transcendental equation. Therefore, the possibility of the 
existence of a power-law tail in the PDF can be verified by this analytical calucula- 
tions. The exponent < s < 1 is also able to be obtained in the same manner as 
mentioned in Sec. [21 It is found that when s ^ \, A ^ and therefore 7 and s 
are connected continuously from s>ltoO<s<l and are continuous functions of 
the interaction parameters a, 6, c, d. A typical example of the comparison between 
the theory and the result of Monte Carlo simulation is illustrated in Fig. ^a). A 
very good agreement is achieved on the right side of PDF. As can be seen from 
Eq. (|3-6|) . there have two cases il) c = d and (II) a + 6 = 1 that s and 7 can be 
derived without evaluating the integral A numerically. It can be shown that these 
special cases have the following interesting features. In the case (I), the parameter 
b is totally eliminated from the transcendental equation despite imposing no restric- 
tions on b. This means that the exchange rate of smaller quantities b is unimportant 
for the system and have no effect on 7 and s. In the case (II), the exponent s is a 
continuous function of three interaction parameters a,c,d and a phase diagram of 
s fixed at d = 0.5 is illustrated in Fig. E^b). It is shown that the diagram of AMI 
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is asymmetrically distorted compared with that of SMI.^^ In the hatched region in 
Fig. n^b), a power-law tail disappears because of the divergence of the solution of 
transcendental equation. It is found that the region is unbounded in AMI, while it 
is bounded in SMI.i°) 




Fig. 1. (a) Double logarithmic plot of versus ^ at a = 0.05, h = 0.9, c = 1.5, d — 0.5. Number 

of particles in a Monte Carlo simulation is A'' = 10^ and number of interactions is T = 50A'^. (b) 
Phase diagram of the exponent a in the case a + 6 = 1 at d = 0.5. Solid curves denote contours 
of s. Power-law tails disappear in the hatched region surrounded by dashed curves. 



§4. Effects of randomness in interaction 

Next, we investigate many-body stochastic processes with multiplicative inter- 
actions with randomly varying interaction parameters.^) For simplicity, the form of 
interaction rule is chosen to be the one of SMI here, that is x' = ax+jSy, y' = [3x+ay. 
The difference from SMI is that at each timestep, the interaction parameters a, (3 are 
randomly determined by a given probability distribution p{a, (5). As in the previous 
section, this model is also considered with following the procedures in Sec. HI A 
master equation for PDF is given by 



dt 



OO fOO fOO 



^ fi^it) = I da I d(3 p{a,(3) I dx I dy f{x,t)f{y,t) 
-'0 Jo Jo 



X [5{z - (ax + py)) - 6{z - x)] . (4-1) 



Using the same scaling relations p-2j) . we have the scaled moment inequalities of 
this processes. 
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where 7 = a + /3 — 1 and the overhnes denote the averages with respect to p{a, (3). 
These inequahties mean that the moment pp diverges at the point where the denom- 
inators of upper and lower bounds of fj,p are equal to zero. This indicates that the 
PDF obeys a power-law tail ^{^) ~ l/^^"*"** whose exponent s > 1 is determined by 
the following transcendental equation. 

S'y = '^ + ^-l. (4-3) 

The exponent < s < 1 can be derived in the same manner as mentioned in Sec. ^ 
Interesting results are found by comparing the growth rate 7 and power-law exponent 
s of this model with those of SMI whose interaction parameters are given by a and 
p. The transcendental equation of this SMI is described by 

sY = a"' +^'' - 1, (4-4) 

where 7' and s' are the growth rate and power-law exponent of the SMI. From 
Jensen's inequality, it immediately follows that a** > a* (s > 1) and a** <a^ {0 < 
s < 1). This concludes that 

7 < 7' (0 < s < 1), s < s. (4-5) 

These inequalities indicate that randomness in interaction parameters have effects 
of reducing the values of both the growth rate and power-law exponent in general. 
In the next subsection, we give one concrete example to check the effects. 

4.1. Two-peaked system 

In this example, p{a, (3) is given by 

p(a, f3) = p5{a - ai)6{(3 - (3i) + (1 - p)6{a - a2)Sif3 - h), (4-6) 

where < p < 1 and ai, /32 are certain positive constants. Here, for conve- 
nience of explanation, we define an index representing the displacement from SMI 
A as follows. 

A = {1- p){a2 - ai)/ao = (1 - p){p2 - Pi)/Po, (4-7) 

where oq = pa\ + (1 — p)a2) /3o = vi^\ + (1 ~ V)f^i- The magnitude of A indicates 
the strength of randomness. The theoretical curves of growth rates and power-law 
exponents are illustrated in Figs. [21 It is obvious that all of their values decrease 
monotonically as A becomes large. It is shown in Fig. Ufa) that the exponents 
drastically decrease even though A slightly increases. Surprisingly, it can be seen 
in Fig. I2tb) that the growth rates reduce from positve to negative when A is large 
enough. 

§5. Summary and discussions 

In this paper, we reported two extended models of many-body stochastic pro- 
cesses with multiplicative interactions from the standpoints of effects of asymmetry 
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Fig. 2. Theoretically evaluated curves of (a) power-law exponents a and (b) growth rates 7 in 
two-peaked system at p = 0.9 (solid lines), p = 0.5 (dashed lines), p = 0.1 (dash-dotted lines) 
with ao = 1.01 and 60 = 0.01. 

and randomness in interaction. It was found by investigating moment inequalities 
that their PDFs have power laws at the tails in general. We successfully obtained 
the values of growth rates and power-law exponents in both models. 

In econophysics, several reserchers^^'^^'^^^'^^^ think that these processes are one 
of the possible candidates for explaining power-law tails observed in wealth distribu- 
tions which arc often called Parcto law^^^ or Zipf's law.-"^^) In economic phenomena, 
there is a well-known concept called the Matthew effect, which is also said "the rich 
gets richer and the poor gets poorer." This effect can be interpreted as a potential 
asymmerty in wealth exchange interactions between transacters. Therefore, it is 
more natural to use AMI for explaining wealth distributions. Remember that in a 
special case of AMI, the growth rate and power-law exponent do not depend on the 
interaction parameter of small quantities. This result would imply that the shapes 
of wealth distributions can be determined by rich people and can be unaffected by 
poor people. On the other hand, everyone has own strategy to transact its wealth in 
real situations. This can be the origin of randomness in interaction. Therefore, it is 
also more natural to consider effects of interaction randomness in the processes. As 
can be seen from the example of two-peaked system, only a few amount of particles 
that has large quantities is sufficient to decrease the values of the growth rate and 
power-law exponent enormously. This result would imply that the entry of a few very 
rich people into a economic system ends up with drastically reducing the economic 
growth rate and magnifying the degree of wealth inequality of the system. These 
implications seem to give us new perspectives on the analysis of wealth distributions. 
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